Dispersion relations and subtractions in hard exclusive processes 
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We study analytical properties of the hard exclusive process amplitudes. We found that QCD factor- 
ization for deeply virtual Compton scattering and hard exclusive vector meson production results in 
the subtracted dispersion relation with the subtraction constant determined by the Polyakov- Weiss 
D-term. The relation of this constant to the fixed pole contribution found by Brodsky, Close and 
Gunion and defined by parton distributions is studied and proved for momentum transfers exceeding 
the typical hadronic scale. The continuation to the real photons limit is considered and the numer- 
ical correspondence between lattice simulations of D-term and low energy Thomson amplitude is 
found. For sufficiently large t the subtraction may be expressed in the form similar to suggested 
earlier for real Compton scattering. 



I. INTRODUCTION 

Hard exclusive reactions described by the General- 
ized Parton Distributions (GPDs) 0, 0, 0, i, 0, H are 
the subject of extensive theoretical and experimental 
studies. The analytical properties of deeply virtual 
Compton scattering (DVCS) and hard exclusive vec- 
tor meson production (VMP) amplitudes @, H, @, [Tol | 
constitute the important aspect of these studies. They 
play also the major role in the hadronic processes such 
as the nucleon-nucleon scattering at very high energies 
to be studied at LHC and in non-accelerator experi- 
ments 



11] 



The crucial point in application of the relevant dis- 
persion relations is a possible ambiguity due to the 
subtraction constants which are the counterparts of 
the normalization constants implied by the ultraviolet 
renormalization procedure. An attractive possibility 
is represented by cases when such constants arc de- 
fined by the imaginary part of the amplitudes. This 
situation was explored long ago in the case of forward 
Compton amplitude 13, 13 . 14 1, and was recently re- 
considered for DVCS [Hit 

In this paper we address the problem of disper- 
sion relations and subtractions in the framework of 
the leading order QCD factorization. We find that it 
leads to subtracted dispersion relations with the sub- 
traction constant defined by the Polyakov- Weiss D- 
term [15j . At the same time, for t exceeding typical 
hadronic scale we relate the subtraction constant to 
the integrals of parton distribution at zero skewness. 



II. DISPERSION RELATION IN THE 
SKEWNESS PLANE 



We restrict our study by the case of large s and 
Q 2 and small t -C s, Q 2 , where QCD factorization 
is applicable. At the leading order, this results in the 
following expressions for DVCS and vector (p°) meson 



production amplitudes: 
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with the GPDs part (which may be interpreted as a 
weighted handbag diagram, i.e. the coupling of local 
quark currents to two photons) 



Af(£, t) 



dx 
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and the meson part for VMP case 
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In 0, H i+ \x, f, t) denotes the singlet (C = +1) 
combination of GPDs, summing the contributions of 
quarks and anti-quarks and of s- and it-channels: 



H f + \x, £, t) 



H f (x,Z t t)-H f (-x,Z,t). (5) 



For the sake of brevity, we will keep only the depen- 
dence of GPDs on the skewness £ = Q 2 /(2s + Q 2 ). For 
£ — * 0, ([3} may acquire divergencies at x = 0, which 
will be one of objects of our analysis. 

Contrary to the forward case, expression does 
not have a form of the dispersion relation because of 
the appearance of £ in the numerator. Nevertheless, 
the amplitude ([3]) as a function of £ manifests the an- 
alyticity in the unphysical region |£| > 1 Q. This 
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region is associated with the contribution of General- 
ized Distribution Amplitudes (GDAs) [l6| related to 
GPDs by crossing [13] . To prove the analyticity of the 
amplitude for |£| > 1 , one represents the denominator 
of (J3j) as the geometric series: 



A(0 
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This series is convergent^ |__thanks to the polynomi 
ality condition (see e.g. 



j dxx n H(x, £) = ^ A 

•J 1, n o 



l-(-l)%n+l 



£ n+i c. 
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One may now easily calculate the discontinuity 
across the cut — 1 < £ < 1 and write the fixed-i dis- 
persion relation Q for the leading order amplitude ([3J 
in the skewness plane: 



V f , DisaA(x) A _, 
RcA(0 = — J dx x _ y ^ + A(£). 
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or, using 
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z-£ 



= V dx 



fl % a) 



A(0, 
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where A(£) is a possible subtraction. This expression 
represents the holographic property of GPD: the rel- 
evant information about hard exclusive amplitudes in 
the considered leading approximation is contained in 
the one-dimensional sections x = ±£ of the two di- 
mensional space of x and £. These holographic as well 
as tomographic (l7| properties in momentum space 
are complementary to the often discussed holography 
and tomography in coordinate space Q . 

We are now going to prove that A(£) is finite and 
independent of £, i.e. A(£) = const. To do this, one 
considers the following representation: 



m =Vld X H(+) ^^-f +) ^^= (9) 



subtraction constant. This can also be derived with 
a use of the Double Distributions (DDs) formalism. 
Namely, the -ff( + )(x, £) is expressed through the cor- 
responding DDs as 



H M (x, 0= J d 



i l-M 

j df3 

i -l+M 



f(a,0)+£g(a,P) 



6(x-a~ £/3) - 5(-x - a - £/3) 



(10) 



Substituting this expression into © , one gets that the 
f(a, /3)-terms which depend on £ are cancelled and 
cqn. J9]) becomes ^-independent (cf. @): 
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A(£) 
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= A. 
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Notice that the cancellation of f(a,/3) and validity 
of (fTTj) are not spoiled (provided £ ^ 0), even if the 
singularity corresponding to f(a,0) ~ a~ a is present. 
In (TTTj) . one can choose 



g(a,0) = 6(a)D(p), 



(12) 



where the function D(/3) is D-term [15|. The assumj. 
tion (fT2"jl is a result of the corresponding "gauge" 
as discussed also in Q. With (fT2"|) . the A term takes 
the following form: 



d/3 



0-1' 
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It should be emphasized that both integrals in © 
are divergent at £ = t = 0, and this divergencies do 
not cancel for £ — > 0. It means that A is not defined 
for £ = t = 26] . On the other hand, for an arbitrar- 
ily small £ the integrals in |(HJ) are finite and, therefore, 
A is well-defined. Note some similarity between (TIU]) 
and ([3]) so that A may be also interpreted as a con- 
tribution of local two-photon coupling to the quark 
currents which is independent of £. 

Taking into account the parameterization [l5| 

oo 

D(0) = (l-0^)Y / dnC^l(0), (14) 



-V [ dxY-^Ht+)(x, £) 
J 4^ n\ d£" 



n-l 



and keeping only the lowest term, one gets 

A = -Ada. (15) 



Due to the polynomiality condition the only surviving 
highest power term in this series is equal to a finite 



This lowest term do was estimated within the frame- 
work of different models. We focus on the results of 
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C M (i\r,) 



chiral quark-soliton model [TP - ^ C Q M ' 
dp = — jjj, where Nf is the number of active flavours, 

and lattice simulations dg att = 
The subtraction constant varies as 



a 



= -0.5. 
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for the DVCS on both the proton and neutron targets. 



III. DISPERSION RELATION IN THE v 
PLANE 



Now, using © with £ = 0, one gets: 



f^^Mi^M (21) 



, H(x, 0) - £/Tx, x) 
IV I dx K ' - J — L , (22) 



where the symmetry property arising from the T- 
invariance: H (x, — x) = H (x, x) is used. The relation 
?T1) can also be obtained from the "sum rules" UJl: 



We have seen that the D-term determines the finite 
subtraction in the dispersion relation in the skewness 
plane. Let us now compare the dispersion relation |8]) 
with the dispersion relation written in the v plane [8j 
where v = (s — u)/4totv. In terms of the new variables 
v' , v related to x, £ as 



x 



4m atz/ 



r 1 



(17) 



the fixed-i dispersion relation becomes the subtracted 
one: 



ReA(v, Q 2 
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Here, = Q 2 /4mjv (the Q 2 dependence here is 
shown explicitly) and the nucleon pole term residing 
in this point may be considered separately Q . 

This subtracted (in the symmetric unphysical point 
v = 0) dispersion relation is the principal result of our 
paper. It i s ap plicable for both DVCS (cf. 0,1) and 
VMP (cf. [if) amplitudes. 

It can be considerably simplified provided lu\A{v) 
decreases fast enough so that both terms in the 
squared brackets can be integrated separately: 



RcAM = — / dv' 2 — 
it J v 



(19) 



where 
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dx 



H(jc,Z + xz)-H(x,£) 



oo ,. 

5^ z n I dxx n - 1 D(x) 
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for £ = and z = 1. 

Let us stress that for the valence (C = —1) contri- 
butions to the amplitudes of the hard exclusive pro- 
duction of, say, pions Q and exotic hybrid mesons 
[2l| . A = because of the mentioned symmetry in x 
and the ^-independence. 

Substituting (|2"Tj) into (|2"0|) , one can see that the 
D-term is cancelled from the expression for the sub- 
traction constant 



C (f) = 2V / dx 



H(x, 0, t) 



(23) 



where we restored the dependence on t which was 
omitted for brevity. This constant is similar to the 
result obtained in the studies of the fixed pole contri- 
bution to the forward Compton amplitude [121 ]. At 
the same time, at £, t = 0, GPDs are expressed 
in terms of standard parton distributions H(x, 0) = 
q(x)6(x) — q(— x)8(— x). Formally one has 



C (0) = 2 dx 



q{x) + q(x) 



dx q v (x) + 2q(x) 



However, the integral defining C o (0) diverges at low 
x in both the valence and sea quark contributions. 
Therefore, for t = we should consider (fT8|) as a cor- 
rect general form of the dispersion relation which in- 
cludes the infinite subtraction at the point v = and 
the subtraction constant associated with the D-term. 

For t ^ 0, the integral in (f23|) converges for suffi- 
ciently large t. In the case of Regge inspired parame- 
terization [4j H(x, 0, —t) ~ x -a{o)+a t ^ £ n j g conc jition 
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reads as t > a(Q)/a', resulting in t > l(10)GeV 2 for 
the valence (sea) quark distributions. 

The divergence of (f2"3"| was originally discussed 
within the framework of the parton model and its 
modifications [l2[ while we address this problem in 
the framework of the leading order QCD factorization. 
Therefore, our result (fT8|) . although being formally £ 
independent, cannot, generally speaking, be contin- 
ued to the forward limit £ ~ Q 2 — > (s = const). The 
limit £ — > 0, s — > oo, Q 2 = const still corresponds to 
the highly non-forward kinematics with the masses of 
initial and final photons being rather different. At the 
same time, further exploration of the possible mani- 
festation of the D-term in the forward Compton scat- 
tering seems very interesting. 

The formal continuation of l|18p would result in the 
following subtracted dispersion relation for the for- 
ward Compton scattering amplitude: 



v 



Rc^(^) = — V 

7T J V 




dv' 2 ImA(v') 



{v' 2 -v 2 ) 



A. (25) 



Comparing this expression with the dispersion rela- 
tion for the forward Compton scattering amplitude 
[22| . one can observe an interesting numerical co- 
incidence. For the proton target, our subtraction 
combined with the lattice simulations (|T6|) is rather 
close to the low energy Thomson term (note that 
Axhomson = 1 f° r °ur normalization of the Compton 
amplitude). As a result, the mysterious occurrence of 
Thomson term at large energies [22j can now be sup- 
plemented with its possible appearance also at large 
Q 2 . This does not hold in the case of neutron tar- 
get where DVCS subtraction term is the same, while 
Thomson term is zero. This may be because the sum 
of squares of valence quark charges is equal to the 
square of proton charge, i.e. to the square of their 
sum (cf. [23l| ) . whereas for the neutron these quanti- 
ties differ. 



IV. CONCLUSIONS 

In this paper we show that the fixed-i dispersion 
relations for the DVCS and VMP amplitudes require 
the infinite subtractions at the unphysical point v = 
with the subtraction constants associated with the D- 
terms. However, for the productions of the mesons 
defined by valence (C = — 1) GPDs the finite subtrac- 
tion is absent. 

We also show that the appearance of the subtrac- 
tion expressed in terms of (forward) parton distribu- 
tions 12| may be investigated in the framework of 
the leading order QCD factorization. We consider the 



possibility of continuation of our results to the real 
photons limit. The surprising similarity between the 
lattice simulations of the D-term and the low energy 
Thomson amplitude in the proton target case is found. 
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